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Abstract 

We present here a study of the radial modes for infinitely long cylinders and cylindrical shells 
of general anisotropic crystals. The elastic coefficients matrix includes twenty-one independent 
constants. We obtain expressions in closed form for the dispersion relation, valid for any anisotropic 
material. In the case of the lowest breathing mode of a thin cylindrical shell we obtain a simple 
analytical formula. This can be used to obtain a first estimate of the breathing mode frequency in 
nanotubes for any material. 
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I. INTRODUCTION 



Cylindrical systems are frequently used as structural components in many engineering 
areas: aerospace, civil engineering, etc., and their vibration characteristics are obviously 
important for practical design. 

The study of acoustic wave propagation in infinitely long homogeneous cylinders and 
cylindrical layers is a well established field from the initial studies of Pochhammer- and 
Chree 2 ' 3 . They developed exact solutions for torsional and longitudinal vibrations of in- 
finitely long homogeneous isotropic solid cylinders with stress free faces. Gazis provided 
solutions for the acoustic wave propagation problem in homogeneous isotropic cylinders of 
infinite length^. Studies for orthotropic hollow and thick cylinders of infinite length, includ- 
ing nine independent elastic constants, were done by Mirsky 5 . Detailed discussions on these 
problems and more references can be found in Refs.6-8. 

No general expressions for the dispersion relation of the acoustic waves in cylinders of 
general anisotropic crystals are available, to our knowledge, although it has been possible to 
obtain solutions in closed form for some particular kinds of modes and crystal systems. The 
radial (breathing) modes of cylinders have a much simpler form than a general vibrational 
mode, thus making them good candidates to try to obtain its dispersion relation for general 
anisotropic crystals. 

On the other hand the study of carbon nanotubes (CNT) and nanowires of many other 
materials, belonging to different crystal systems, has increased the use of isotropic con- 
tinuum models for cylinders and cylindrical shells^ - —. In the case of transverse elastic 
isotropy (hexagonal crystals) formal expressions were given in Ref.14. In the case of general 
anisotropic crystals, a method developed in resonant ultrasound spectroscopy to get the free 
vibrational modes of inhomogeneous systems^— , has been employed in several studies of 
acoustic modes in nanowires^"—. 

Quite recently the acoustic modes in nanorods of different materials 2 ^— have been ob- 
served directly by optical methods. Other approaches in continuum medial— are also 
been used. This shows the interest in the study of the vibrational modes of nanotubes and 
nanowires. 

Thus it is worthwile to try to get expressions in closed form for the radial modes frequen- 
cies in cylinders and cylindrical shells of general anisotropic crystals. 
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We shall consider here an infinitely long cylinder, or cylindrical shell, of a general 
anisotropic material and we shall obtain the dispersion relation of the radial modes for these 
systems. 

In Section II we present the formal equations leading to the dispersion relation and the 
comparison with other simpler cases. Conclusions are presented in Section III. 

II. ACOUSTIC RADIAL MODES IN CYLINDERS AND CYLINDRICAL SHELLS 
OF GENERAL ANISOTROPIC CRYSTALS 

We shall consider a general anisotropic crystal. The matrix of elastic coefficients will by 
given in this case by 
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C14 


C*15 


C*16 




C*22 


C*23 


C*24 


C*25 


C*26 


C*13 


C*23 


C33 


C34 


C35 


C*36 


C14 


C24 


C34 


C44 


C45 


C*46 


C*15 


C*25 


C35 


C45 


C55 


C*56 


Cl6 


C26 


C36 


C46 


C*56 


C*66 



(1) 



including twenty-one independent constants. 

We shall consider a homogeneous thick cylinder of infinite length and radius R, or a 
homogeneous cylindrical shell of infinite length, inner radius a, outer radius b and thickness 
h (b = a + h). We shall use the cylindrical coordinate system (r,9,z). 

We assume that the mass density p is constant in the cylinder (cylindrical shell). The 
cylinder axis coincides with the crystalline axis z and the surfaces are stress free. 

The strain-displacement relations are 



du r 1 / dua \ du 



or r \ 06 J Oz 

1/1 du r dug Ug\ 1 / dug du z \ 1 / du r du z ■ 

er9= 2\r^e + ^r~~VJ ' ,€ez= 2 \8z~ + ~d6 ) ] " rz = 2 V ~dz~ + ~dr~ ' ' (2) 

The equations of motion are given now by 
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(3) 



where 



a = C e . (4) 

The general solution of this problem, as we told before, must be performed numerically, 
but we shall consider now the particular case of the radial modes. 

The radial modes are purely radial vibrations, so that u r ^ and u e = u z = 0. In the 
same way 

du r du r 

It can be shown that the second and third equations in (3) are satisfied automatically 
and we are left with 

d " Ur 1 r^r + (A 2 - = o , (6) 



dr 2 r dr 1 r 

2 _ ^22 



r, OOO 2 9 C' 

with pf = — — and /T = — . 
Eq.(6) is a Bessel equation of non-integer order. Thus we have the solution 



AJ^ftr) + BJ^ifrr) , ft > 
u r (r) = ' B (7) 

Ar + - , A = . 

r 

The boundary conditions at the surfaces are now 

C n ^ + C 12 - = 0. (8) 

dr r 

In the case of a thick cylinder, in order that u r be finite at r=0, B—0. Thus we shall 
have: 
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AJ^fyr) , A > 
u r (r) = ^ (9) 



If cr rr (i?)=0, A=0 when A=0, and when A > we have 



As 



CnJ^R) + Cu-^MftR) = . (10) 



j;(Ar) = A^-i(Ar) - ~Jn(Pir) (11) 



eq.(10) becomes 

1 

7? 

thus giving 



CnPiJ,*-i(PiR) + (Ci 2 - Cnri-J^ftR) = , (12) 



ftRJ^ftR) = (fi- ^)UfrR) ■ (13) 



This equation is a generalization of the expression obtained for an infinitely long isotropic 
cylinder 29 

ftMftR) = ^\~ 2U) J^R) (14) 

Q 

with the Poisson's factor v 



2(Ci2 + C44) 

In the isotropic case Cn = C 2 2 and Cn = Ci 2 + 2C44. Thus it is easy to see that eq.(13) 
becomes in that case eq.(14). 

In the case of the cylindrical shell, by using the boundary condition given in eq.(8), we 
obtain: 



A(C u j;(Ar) + cJ-J^r)) + S(C u J^(Ar) + C 12 l - J-^&r)), A > 

A{C lt - C X2 ) - B{C 1X - C 12 )i A = . (15) 

As a rr (a) = a rr (b)—0, we must have A = B=0, for A=0, whereas for A >0 we have 



CnA^(Aa) + Cu-Mfca) C^JL^ia) + C 12 - J-^fta) 



CnA^(A6) + Cia^Afe) C^J'^b) + C 12 ^ J_ M (A&) 



= . 



(16) 



This is the dispersion relation for the elastic radial modes of an infinitely long cylindrical 
shell in the case of general anisotropy. In hexagonal crystals C 22 = Cn and C\i = Cn — 2C 66 . 
Then the last equation reduces to eq.(33) in Ref.14, for the hexagonal crystals. 

In the case of a thin shell it is possible to make an expansion to first order in — in eq.(16). 

a 

Then we arrive to 



[(^ 2 -CnC 22 )+Cj^ 



(17) 



Taking into account that 



a 



(18) 



it can be seen that the factor involving the Bessel functions in eq.(17) is the Wronskian, 
having the value 



U(3 ia )J'_,((3ia) - J^(Aa) J_„(Aa) = _£E^?0 ^ . 



(19) 



Thus the dispersion relation for the lowest breathing mode of a thin cylindrical shell of a 
general anisotropic crystal is given by 



(C? 2 - C n C 22 ) + C^flV = , 



(20) 



which can be put in the form 



1 / CllC22 — C\2 



2 / CnC 22 — C\ 2 



pC n 



d 



(21) 



d = 2a being the diameter. 
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In the case of hexagonal crystals C22 = Cu and C 12 = Cu — 2C 66 and eq.(21) reduces to 



which coincides with the expressions given in Refs. 11,14. 

Eq.(21) is a simple analytic expression valid for all crystal systems. It provides in a quick 
way a first estimate of the lowest breathing mode frequency for nanowires. It was shown in 
Ref.14 that eq.(22) gave extremely good agreement with experimental and first principles 
theoretical values of the lowest breathing mode frequency of different nanotubes. It can be 
expected that the same will happen for eq.(21) when applied to other materials and crystal 
systems. 

Eq.(16) gives the frequencies of higher order breathing modes. 

III. CONCLUSIONS 

We have obtained expressions in closed form for the dispersion relation of radial modes 
in infinitely long cylinders and cylindrical shells of general anisotropic materials. When 
considering the case of a very thin shell we arrive to a very simple analytic expression for 
the lowest breathing mode. This expression covers the cases of isotropic and hexagonal 
systems, previously given in the literature. 
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